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Abstract .  'Acceleration  waves  in  elastic-plastic  materials  are  studied  in 
some  detail  on  the  basis  of  a  nonlinear  thermodynamical  theory  of  elastic- 
plastic  continua.  Attention  is  confined  mainly  to  non-conducting  media, 
but  the  developments  are,  otherwise,  general.  Formulae  for  wave  speeds 
are  derived,  fronts  of  plastic  loading  and  elastic  unloading  are  discusse 
and  higher  order  discontinuities  are  shown  to  have  the  same  characteristic 
speeds  as  those  of  acceleration  waves .  An  example  concerning  propagation 
of  plastic  waves  in  a  medium  undergoing  uni-axial  motion  is  included. 
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Is  an  independent  thermodynamic  variable  and  then  obtain  an  expression 
for  the  rate  of  production  of  entropy  in  a  non-conducting  medium*  It  is 
a  known  result  in  nonlinear  elasticity  that  in  a  non-conducting  elastic 
material  (and  in  the  absence  of  heat  supply)  the  time  rate  of  entropy 
vanishes  and  the  so-called  acoustic  tensor  is  homentropic.  In  general, 
as  shown  in  sections  3  and  4,  such  is  not  the  case  fo v  an  elastic-plastic 
material. 

In  section  4,  waves  are  represented  as  propagating  surfaces  in  a 
non-conducting  elastic-plastic  continuum,  across  whicn  there  exiso  jumps 
in  the  values  of  certain  mechanical  and  thermal  variables.  Acceleration 
waves,  or  singular  surfaces  with  respect  to  particle  acceleration,  are 
discussed  in  detail  and  formulae  for  wave  speeds  are  derived  using  the 
compatibility  conditions  of  Kaaamard  [13^  for  jumps  in  piecewise  continuous 
functions  of  position  and  time .  Waves  propagating  through  regions  under¬ 
going  plastic  deformation,  sufficient  conditions  lor  the  existence  of  re? a. 
wave  speeds ,  and  fronts  of  plastic  loading  and  elastic  unloading  are 
treated  in  sections  4-6.  Singular  surfaces  with  respect  to  any  order 
derivatives  (in  time,  position  or  mixed)  of  acceleration  are  shown  in 
section  7  to  have  the  same  characteristic  speeds  as  tnose  oi  the  accelera¬ 
tion  discontinuities.  Finally,  in  section  8,  we  consider  the  simple 
example  of  propagation  of  plastic  waves  in  a  medium  undergoing  uni-axial 
motion.  Other  examples,  such  as  plastic  waves  in  simple  shear  or  propa¬ 
gation  of  spherically  symmetric  plastic  waves,  can  be  discussed  in  a 
similar  fashion. 

The  above  developments  for  a  non-conducting  medium  in  sections  3~7 
are  valid  for  a  work-hardening  elastic-plastic  material  which  is  initially 
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anisotropic  and  are  ad. so  applicable  to  the  limiting  case  of  an  elastic- 
perfect.ly  plastic  material.  Of  special  interest  is  the  simple  formula 
that  the  jump  in  the  entropy  production  is  linear  in  the  jump  of  plastic 
strain  rate.  This  and  other  results  in  section  4  hold  also  for  a  definite 
conductor,  i.e.,  for  a  medium  whose  heat  conduction  vector  has  the  form* 
of  Fourier's  law. 
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2 .  Notation  and  Preliminaries 


L°t  the  motion  of  the  continuum  be  referred  to  a  fixed  system  of 
rectangular  Cartesian  axes  and  denote  the  position  of  a  typical  particle 
at  time  t  by+ 


x. 

i 


(2.1) 


where  X ^  is  a  reference  position  of  the  particle.  Vie  require  the  mapping 
(2.1)  to  be  single-valued  and  have  continuous  partial  derivatives  with 
respect  to  its  arguments,  except  at  some  singular  points,  curves  and 
surfaces.  We  use  the  notation  F  =  F(t)  and  designate  partial  dif¬ 
ferentiation  with  respect  to  X.  or  x.  as  (  )  or  (  )  .  , 

A  2-  3  A  3  1 

respectively.  Latin  indices  take  values  1,2,3  and,  except  when  noted 
otherwise,  the  usual  summation  convention  for  Cartesian  tensors  will  be 
employed. 

The  components  of  velocity  and  acceleration  at  the  point  x^ 
at  time  t  are,  respectively, 


x. 

i 


h  *i(xA’ 


•  # 
x. 

i 


_d_ 

at 


Xi(x 


A’ 


t) 


(2.2) 


where  a  superposed  dot  stands  for  differentiation  with  respect  to  t  , 

holding  X.  fixed.  We  define  the  strain  tensor  e  by 
A  KXj 


+We  use  the  same  symbol  for  a  function  and  its  value  without  confusion. 
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6. 


heat  flux  vector,  both  zensnrad  per  emit  tics  sac  per-  cr-it  area  in.  the 
reference  configuration.  ahem,  zhs  (lacal)  energy  equation  and  the 
(local)  C3  sasias-lteaas  inequality  are 
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beat  supply  function  per  unit  zrass  and  per  unit  ti^e. 


U  is  the  internal  energy  per  unit  rass,  S  is  tbe  entropy  per  unit 

mass  and  *!(>  0)  is  tbe  temperature* 

The  kinersatical  Quantity  e„  is  unaltered  nben  the  continuum  is 

.u* 

subjected  to  superposed  rigid  body  notions  at  ail  tines  t  =  (tr&)  , 
where  a  is  a  constant.  In  what  follows  we  shall  be  concerned  with 
constitutive  equations  which  mainly  involve  stress,  stress  rate,  dis¬ 
placement  gradient,  velocity  gradient  and  such  equations  must  remain 
unaltered  by  superposed  rigid  body  motions .  As  noted  in  [ll ,  the 
displacement  gradient  and  the  velocity  gradient  must  be  replaced  by 

e„  and  e__  and  these  will  be  referred  to  as  the  strain  tensor  and 
KL  KL 

the  strain  rate  tensor.  Also  the  stress  and  the  stress  rate  may  be 

taken  to  be  s  and  s  ,  both  being  invariant  under  superposed 
KL  KL 

rigid  body  motions. 
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3.  Constitutive  Eauations  and  Some  Related  Results 


We  summarize  here  the  principal  results  from  the  theory  of  elastic- 

plastic  eontinua  as  developed  by  Green  and  Naghdi  [l,2]  and  also  obtain  some 

related,  results  for.  later  use.  We  limit  the  discussion  to  the  material_  de- _ 

script ion  of  the  theory. and,  for  convenience,  regard  the  entropy  S  (rather _ 

than  the  temperature)  as  an  independent  thermodynamic  variable. 

Let  the  strain  tensor  e„.  be  defined  at  each  point  of  the  continuum 

by  (2.3)  and  let  e"  ,  a  symmetric  tensor  with  the  same  invariance 

AL 

properties  as  e.-  ,  denote  plastic  strain.  We  introduce  a  constitutive 

AJLi 

assumption  for  sT,_  in  the  form 
KL 


SKL  “  SXL(eKL’  eKL5  S)  ’ 


(3.1) 


and  admit  the  existence  of  a  scalar-valued  continuously  differentiable 

function  f(s  ,  e"  ,  s)  —  called  a  yield  or  a  loading  function  —  such 
A_Li  Al 

that  the  equation 


f(SKL’  eKl’  "  K  » 


(3.2) 


for  a  fixed  value  of  K  and  e"  represents  a  hypersurface  in  seven- 

IUj 

dimensional  Euclidean  space  --  six  components  of  s„  and  the  entropy 

Ali 

S  .  The  scalar  K  ,  a  work-hardening  parameter,  depends  on  the  past 


’’"‘"Although  we  assume  that  e"^  is  symmetric,  the  developments  which 

follow  can  be  readily  modified  to  accommodate  a  non- symmetric 
plastic  strain  tensor  included  in  [2,143 > 
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postulates  for  the  internal  energy __  U  ,  the  temperature  T  _and  the  i  sat 

* 

flux  vector  Q.^  .  Thus,  we  introduce . 


U  *  U(eKW  ek>  S) 


and  also  assume  that  T  and  Q„  are  functions  of  e„T ,  e"  ,  S  an ,  that 

iv  KL  KL 

a.,  depends  in  addition  on  T  . 

K  ,M 

With  the  foregoing  background,  we  now  record  certain  addition  1 
results  from  the  theory  of  elastic-plastic  continua  [1,2]  whica  w_  1.1  be 
utilized  subsequently.  The  constitutive  equations  for  e1'  ar.d  /.  are 

nii 


'KL 


K  =  h 


3KL 


KL  KL 


A 

both  of  which  hold  during  loading,  i.e.,  when  f  =  K  ,  f  >  0  or  during 

•  A 

neutral  loading,  i.e.,  when  f  =  K  ,  1C  =  0  and  f  =  0  .  In  <3«'i )  and- 

(3.8) ,  >*  and  p  are  tensor  functions  and  i  is  a  positive  scalar 

KL  KL 

function  of  s,_ ,  e"  and  S  .  In  addition,  we  have 
ML  ML 


* 


As  in  [2] ,  we  may  allow  U,  T,  Q„  ,  as  well  as  the  stress  ter  or  s  , 

K  KL 

to  depend  also  on  K  whose  rate  is  specified  by  (3.8)  below  Vor 

simplicity,  we  retain  here  the  constitutive  assumptions  of  t  ■ ;  forms 

(3.6)  and  (3.1)  but  note  that  the  inclusion  of  K  as  an  independent 

variable  will  not  alter  the  structure  of  most  of  the  results  m  later 

sections  of  the  paper. 
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which  are  independent  at  rates  and  are  valid  both  during  loading  and  un- 
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idsrstood  to  have  the  symmetric  form 
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*5 ee  £reen  and  laghdi  [li  for  details.  Although  the  main  developments 
in  [li  are  carried  cut  in  terms  of  the  Helmholtz  free  energy  function 
(with  r  as  an  independent  thermodynamic  variable),  results  cor¬ 
responding  to  (3-6)  and  (3.10)  are  also  discussed  [l;  section  6] .  In 
the  }  inert  paper,  as  in  [2J ,  the  basic  kinematic  variables  are  e^ 

which  is  defined  by  (2.3),  while  the  original  form  of  the 

Ln  [li  enuloyed  ef_  and  a  variable  defined  by  e  '  =  e„  -  e"  . 

"  <Uj  iUi  AJ-)  iiL 


ariz  a*. 


Ihe  two  forms  of  the  nonlinear  theory  in  [1,2]  are  entirely  equivalent, 
although  in  general  the  latter  [2]  is  preferable.  As  emphasized  in 
rll,  the  variable  e'  is  not  an  "elastic"  strain  in  the  usual  context 

*"  7  fCT, 

of  elasticity;  only  in  restrictive  cases  or  in  the  infinitesimal  theory 
of  el asti c-nlasti c  continua,  e '  is  an  elastic  strain  tensor  ana  is 
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where  (3.13)  is  deduced  frora  (3-12)  by  considering  an  arbitrary  homogeneous 

temperature  distribution  for  which  T  ,,  =  0  and  recalling  chat  X  >  0  . 

jk 

The  above  results  are  valid  for  any  elastic-plastic  continuum.  For 


later  reference,  we  note  that  a  medium  is  called  elastic -perfectly  plastic 


if  the  loading  function  f  and  the  tensor  function  reduce  to  the 

-r*f 

iorms 
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Ko  ’  ^KL  ~  ^KL^SMN’  ’  (3*15) 

where  Kq  is  a  real  constant.  Since  f  is  now  independent  of  e^_  , 
the  loading  surface  is  always  stationary.  In  this  case,  all  terms  in- 

— S  .p 

volving  --°v~  vanish  and  h„.  in  (3-8)  must  also  vanish.  Neutral 
3eKL  151 

loading  no  longer  exists  and  the  condition  for  loading  reduces  to 
A 

f  =  0  .  Furthermore,  it  follows  from  (3*9)  that  for  an  elastic-perfectly 

plastic  medium  X  -*  »  with  87rT  remaining  finite  and,  in  place  of  (3-7), 

iLL 

we  have 


(3-16) 
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eia  ^kl 

PPQ  ^PQ, 


(3.17) 


‘  The  definition  of  an  elastic-perfectly  plastic  material  here  is  not 
the  same  as  that  in  [l;  section  9) >  obtained  by  specialization  of  the 
spatial  form  of  the  theory.  A  number  of  different  definitions  of  elastic- 
perfectly  plastic  materials  are  possible  in  the  context  of  the  nonlinear 
uheory,  but  they  coincide  for  a  linearized  theory. 


The  results  (3-12)  and  (3.14)  remain  valid  for  an  elastic-perfeetly 
plastic  medium,  hut  the  inequality  (3.13)  must  he  replaced  by 
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In  this  paper,  we  restrict  our  attention  mainly  to  a  medium  which 

is  a  non-conductor  and  also  exclude  the  heat  supply  arising  from  external 

>  •  ^ 
neat  sources  and  energy  losses  due  to  radiation,  i.e., 


0  ,  r  a  0 


Then,  from  (3.12),  we  get 
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Also,  hy  (2.9), 
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(3-1$) 


(3-20) 


(3.21) 


Remembering  the  constitutive  assumption  (3-6),  we  differentiate  (3.10)^ 
to  obtain 


It  is  not  essential  to  assume  '3.19)2  which  is  introduced  here  for 

simplicity  of  subsequent  formulae.  Further  comment  on  this  is  made 
in  section  4. 
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3_^  is  rcn-zero.  Sert,  substitute  (3-26)  into  (3-23) 

£r  siuiiar  uc  bbat  xhich  led  to  (3-26)  solve  for  s  in 
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(3-26)  and  (3-27). 
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k.  Acceleration  Discontinuities 


-Let  (X, ,  t)  denote  the  position  at  time  t  of  the  material 

1  1  xi 

■scant  whose  position  In  the  reference  configuration  is  X^  ana  consider 

a  moving  surface  in  the  continuum,  a  smoouh  one-parameter  family  of  points, 

across  whi oh  certain  derivatives  of  x,.  (Xfl  >  u)  have  jump  discontinuities  * 

Such  a  surface  is  called  a  singular  surface  or  a  wave  and  may  he  assigned 

*)  2  1^2 

the  material  description  v(t)  :  X,  =  Yfi(8~,  S  >  t)  ,  where  6'  and  e 

,1.2  .. 

are  parameters.  Alternatively,  hy  elimination  oi  8  ana  6  ,  "one 

smooth  surface  £(t)  may  he  represented  as 


S(t)  :  5(Xa,  t)  -  0  .  fr-D 

Ecuatior.  (4.1)  locates  the  surface  as  a  function  of  time  in  the  reference 
configuration . 

The  urdu  normal  to  X(t)  will  he  denoted  hy 
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jne  speed  of  propagation  of  the  surface  (4.l)  is 


(4.3) 
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Let  R  and  RT  be  one  sided  adjacent  neighborhoods  partitioned 
by  s(t)  and  let  F(X^,  t)  ,  a  function  of  X^  ,  be  continuous  in  R- 
ana  R+  but  have  a  jump  across  x(t)  .  We  define  the  jump  in  k(X£ ,  t) 
across  2  at  time  t  by 

[[f|=  F+  -  F~ 


=  F(Y+(e*,t)}t)  -  F(Y“(ec',t),t)  , 


(4.5) 


where  ?'  ana  F  designate  the  values  of  F  on  the  two  sides  of 
2(t)  approached  from  the  neighborhoods  R'  and  R  ,  respectively. 

in  what  follows,  we  mainly  consider  a  singular  surface  of  order  2, 
i.e.,  an  acceleration  wave  but  also  briefly  discuss  singular  surfaces  of 

■jf  .  . 

order  m  >  2  .  We  call  2(t )  an  acceleration  wave  in  elastic-plastic  con- 
■  cinua.it'  xn,  x^  ^  (and  therefore  e^) ,  e"_  and  S  (or  T)  are  continuous  func¬ 
tions  of  X^  and  t,  while  their  derivatives  (with  respect  ~o  X  or  t)  have 

at  most  jump  discontinuities  across  £(t)  but  are  continuous  in  R‘  and 
R  .  We  also  assume  that  the  externally  applied  body  force  is  so  assigned 

that  F.  is  a  continuous  function  of  X.  for  all  t  .  It  then  follows 

l  -  A 

•chat  the  kinematical  conditions  of  compatibility  are  (see,  e.g.,  Truesdell 
and  Toupin  [15  >  section  190l) 
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’A  singular  surface  of  order  0  is  a  stationary  surface  and  a  singular 
surface  of  order  1  is  called  a  shock  wave.  The  terminology  used  here  is 
due  to  Hadamard  [13];  see  also  Truesdell  and  Toupin  [15]  or  Thomas  [l63 . 
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ana  equation  (4.9)  beccces 


-  *3 *i,A I^!+  v<Wab  -  Po -  0  »  P*-11) 

where  (4.6)  and  (4.10)  have  been  nsec.  Also,  frca  (3-20),  we  have 

(4.12) 

Provided  BU/de t-  does  not  vanish,  according  to  (4.12),  the  jump  in 
the  entropy  production  is  linear  in  the  jump  of  plastic  strain  rate. 

It  is  therefore  clear  that  acceleration  waves  (with  V  /  0),  in  a  non¬ 
conducting  elastic-plastic  medium,  are  not  hoaentropicT 

In  the  rest  of  this  section,  we  consider  three  types  of  wave  propa¬ 
gation  as  follows: 

(i)  In  the  absence  of  loading  froo  an  elastic-plastic  state,  on 
either  side  of  the  singular  surface  s(t)  ,  we  have  either  a  state  of 
unloading  or  neutral  loading.  Recalling  (3.29),  we  then  write 


It  is  not  difficult  to  see  that  the  result .(4.12)  holds  even  without 
the  assumption  (3.19)25  it  Is  only  necessary  to  assume  that  r  is  a 

continuous  function  of  X,  for  all  t  . 

A 

+An  acceleration  wave  is  called  homentropic  if 
(Recall  that  if  |  S  jj  =  0  ,  it  follows  from  |  S  |J  =  0  that 
![  S  M  1  =  ^ne  relation  (4.12)  may  be  contrasted  with  the  cor¬ 

responding  result  obtained  by  Coleman  and  Gurtin  [5).  For  the  class 
of  materials  with  memory  discussed  in  [5] ,  they  found  that  every  ac¬ 
celeration  wave  (with  V  /  0)  in  a  non-conductor  is  homentropic. 
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and  with  the  help  of  (2.4)  and  (4.6)  this  becomes 
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Using  the  notation  V  =  ^  for  this  case,  substitution  of  (4.14)  into 

(4.11)  results  in 


(A.  .  ”  ( \  —  ““  S/^-\  )§.  .  .  —  0 
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and  the  associated  eigenvalue  equation  is 


det  fA. . 
1  2-0 


(V^  v  -  —  s,„N)6. .]  =  0 
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(4.16) 


where  we  have  defined 


s(n)  skl  nk  nl  5 


A.  .  =  A. .  =  x.  _  x.  „  K„  K..  A.„._.r 
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Mathematically,  this  case  is  similar  to  that  for  nonlinear  elasticity 
with  playing  the  role  of  the  homentropic  acoustic  tensor  and  with 

|  S  J  =  0  .  Here,  however,  the  tensor  A^.  depends  on  the  current 
state  of  -clastic  strain  e"  ,  as  well  as  on  e  _  and  S  ;  it  can  be 

IU_|  XUj 

identified  with  the  homentropic  acoustic  tensor  of  elasticity  only  when 
the  medium  has  no  prior  history  of  plastic  deformation. 
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j;v  .  wzwwj;  p-Tassroro  •■■ 
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-  ’•  /'<:•: -4^* f%5u'7^|^ 

*  .  •  M^fe? 

undergoing  plastic  deformation,  i.e.,  a  plastic  wave,  and  will  je  further  _j  * 


discussed  in  the  next  section. 

(iii)  Finally  we  use  the  notation  V  =  to  correspoi.  .  to  the 

case  when  the  region  RT  is  in  a  state  of  loading  and  the  region  R 
is  in  a  state  of  unloading  or  neutral  loading,  with  V  =  correspond¬ 

ing  to  the  reverse  of  this  situation.  In  the  first  of  these  <<.sss,  it 
follows  from  (3.27)  and  (3.2 9)  that 


1  SXL  1  Po  AKD®  I  ®MN  1  “  Po  GKL  SM 


"_po  V(u)  AKLMH  NN  Xk,M  \  '  po  C'KL  4®  eMW 


•(4.22) 


--  :A  1 


Hence,  with  the  help  of  (4.11),  we  obtain 


^Aij  ~  <4)  "  pQ  a(N);oi3  j/v0  "  wi  1JM  CM? 


S/„v)6.  .}X.  =  -  — —  G_.  H,„t  e* 


(4.23) 


On  the  other  hand,  when  R+  is  in  a  state  of  unloading  or  neutral  loading 
and  R~  is  in  a  state  of  leading,  by  (3*27)  and  (3.29)  we  nave 


!  SXL  1  -  po  AKD®  I  em  ]j  +  po  GKL  4®  eMN 


=  po^AKLMN  ~  gkl  *W  I  eMN  1  *  po  gkl  4®  6 


® 


p0V(^'/Am®"  GKL  *WNN  *k,M  Xk  +  po  GK  4®  eM  *  (4.24) 


'"in  the  literature,  the  term  "plastic  wave"  is  often  usi  1  to  denote  a 
moving  elastic-plastic  interface  (see,  e.g.,  Hopkins  [  ?]). 
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vhcse  heat  conduction  vector  h  the  form 


~Zcuation  (4.27)  is  in  the  form  of  Fourier's  lav;  and  the  heat  conduction 
tensor  YL-  &  function  of  e^.,  e"r.  and  T,  is  positive  definite. 
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(4.27) 


T 

M 


In  this  case,  it  is  more  convenient  to  choose  T  (rather  than  S)  as  an 
independent  thermodynamic  variable  and  to  use  the  form  of  the  theory  in 
which  the  constitutive  equations  are  expressed  in  terms  of  the  Helmholtz 
free  energy  function 


A  =  U  -  TS  ,  (4.28) 

together  with  appropriate  changes  in  section”  3*  In  view  of  the  con¬ 
stitutive  and  smoothness  assumptions,  it  follows  from  the  integral  form 
of  the  equation  of  balance  of  energy  that  every  acceleration  wave,  in  an 
elastic -plastic  material  subject  to  Fourier's  law  of  heat  conduction,  is 
homothermal,  i.e.. 


Moreover,  we  now  obtain 


(4.29) 


(4.30) 


in  place  of  (4.12).  It  should  now  be  clear  that  results  parallel  to  those 
between  (4.13)  and  (4.26)  can  also  be  deduced  in  this  case,  but  we  do  not 
pursue  the  matter  further. 


T'?or  details,  see  [1,2].  In  this  form  of  the  theory,  A,  s^  and  S,  as 

well  as  f,  h,_  and  B._,  are  functions  of  e._T,  e"  and  T.  Also,  instead 
XL  XL  MN  Mn 

of  (3.19)9,  we  only  assume  that  the  heat  supply  is  so  assigned  that  r  is 

a  continuous  function  of  X.  for  all  t . 

This  conclusion  parallel^  the  corresponding  result  for  elastic  materials 
[3]  and  for  the  class  of  materials  with  memory  considered  in  [5]  with 
the  heat  flux  vector  in  the  form  (4.27).  For  a  definite  conductor,  the 
heat  conduction  tensor  in  [5]  is  given  a  more  general  definition  than 
that  in  (4.27). 
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5.  Plastic  vr&ves 


r(P> 

must  identically  satisfy  the  eigenvalue  equation  (4.21)  fear  given  values 


,  a 


'XL5  XL 


and  S  .  Of  course,  one  converse  is  net  true.  i.e..  every 


solution  7^  to  (4.21),  for  give 


an  values  of  e._ .  ef_ .  s  -  does  not 


necessarily  represent  the  speed  of  an  acceleration  wave.  Bevever,  tie 
real  values  of  Y^  are  the  only  speeds  at  vhieh  such  -saves  nay  propagate. 

To  investigate  sufficient  conditions  for  real  wave  speeds,  ve  appeal 
to  known  properties  of  matrices .  Consider  A. .  .  defined  'ey  vheri 

is  given  by  (3*25).  I?  is  an  arbitrary  vector,  then 


£  -s  t>  =  fp  Vp  y  Vlf  v  a 

"ij  -  i  - j  K~  1  i,Ln* 3  j  ,~iJ  £  'V.  "XlMZs 


which,  in  view  of  the  symmetry  properties  of  , 


2.7  ce 


ii  P  P  =  £  S  W 

"ij  ri  j  "XI2‘5f  XL  XS  5 


/r 


VJ_  being  a  symmetric  tensor  defined  by 
aL 


W._  =  ?.  x.  _  3L  x.  „) 

XL  2  iv  X  i,i.  l  l,X 


(5-2) 


From  (5  .1) ,  one  concludes  that  A.  is  a  positive  de: 


is  positive  definite  in  trie  sense  tnat 


A  w  w  a  0 
rCLMN  XL  *JS  “ 


Ky-sJ 
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for  all  arbitrary  symmetric  ,  with  the  equality  holding  if  and  only 

if  is  identically  zero. 

If  (5 -3)  holds,  then  the  matrix  has  three  real  eigenvalues, 

which  may  be  ordered  numerically  as 

5  Ag  Bs  A^  >  (5*4) 

and  which  Identically  satisfy  the  equation 

det  \k±i  -  ka  6ij |  =  C  (a-  1,2,3)  •  (5-5) 

Associated  with  these  eigenvalues  are  a  set  of  three  eigenvectors 
0^c^  (a  =1,2,3)  which  are  found  to  within  a  scalar  magnitude  by  the 
homogeneous  system  of  equations 

CAU  '  Aa  ‘id1  0dO)  '  0  ’  (5-6) 

Now,  following  the  technique  used  by  Mandel  [10],  let  the  eigenvectors 
0^  be  referred  to  the  principal  directions  of  the  matrix  A^  and 
write 


Ai  6id 


(5.7) 


where  a  bar  below  an  index  signifies  suspension  of  the  summation  convention 
for  that  index.  Then,  the  eigenvalue  equation  (4.19)  for  plastic  waves 
becomes 
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(5.3) 


4«t  KAj  -  Y)6lj(  -  Ot  Hjl  .  0  , 

which  when  expanded  hae  the  form 

P(Y)  ■  (A1  -  Y)(Ag  -  Y)(A3  -  Y)  -  G1  H^Ag  -  Y)(A3  -  Y) 

-  02  h2(a3  -  Y)(AX  -  Y)  -  o3  h3(a1  -  Y)(A2  -  Y)  -  0  ,  (5.9) 

where 


Y  ■  4)  ■  i  9(k) 


(5-10) 


Since  (5-9)  is  &  cubic  polynomial,  at  least  one  of  the  roots  Y  (a  ■  1,2, 3) 
must  be  real.  Of  course,  with  the  help  of  (5 .10),  it  is  seen  that  the 
wave  speed  corresponding  to  Y^  will  be  real  only  if  ' 


Y  + 
a 


s  0 


(5.11) 


We  also  find  from  (5*9)  that 


P(A1)  -  -  G1  H1(A1  -  A2)(A1  -  A3)  , 


P(Ag)  -  G2  H2(A1  -  A2)(A2  -  A3) 


(5.12) 


p(a3)  -  -  o3  h3(a1  -  a3)(a2  -  a3) 
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In  view  of  (5.4),  w©  observe  from  (5-12)  that  P^),  P(Ag)  and  P(A^) 
have  the  signs  of  -  G^  ,  +  Gg  Hg  and  -  G^  ,  respectively.  Noting 
that  P(Y)  -  os  as  Y  -*  +  «  while  P(Y)  -*  +  »  as  Y  -*-<»,  we  deduce 
the  following  information  about  the  wave  speeds: 

If  H^,  Gg  Hg  and  G^  all  have  the  same  sign,  then  Y^,  Yg 
and  Y^  are  real  and  distinct;  if  this  sign  is  positive,  then 

Ax  S  Yx  ^  Ag  5  Yg  5  A3  5  Y3  ,  (5-13) 

while  if  it  is  negative, 


I15A1?I2S  A2  *  Y3  *  A3 


(5.14) 


(a) 

It  is  clear,  however,  that  if  is  to  be  a  real  wave  speed,  the 

inequality  (5.11)  must  hold. 

By  comparing  (5*6)  with  (4.15),  we  further  observe  that  A  are 

CY 

precisely  the  eigenvalues  for  the  case  where  there  is  no  loading  (i.e., 
additional  plastic  deformation)  on  either  side  of  the  wave  front .  It 
would,  therefore,  be  of  physical  interest  to  compare  A^  ,  the  speed  of 
the  fastest  "elastic"  wave,  with  the  real  roots  among  Y^,  Yg  and  Y3  , 
Similar  comparisons  for  Ag  and  may  be  of  interest. 
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6-  Fronts  of  Loading  and  Unloading 


Any  surface  in  the  medium  which  separates  a  region  of  loading  from 
one  of  unloading  or  neutral  loading  is  called  an  elastic-plastic  interface. 
With  reference  to  case  (iii)  of  section  4 ,  suppose  such  a  surface  is  moving 
through  the  medium  with  an  acceleration  discontinuity  coinciding  with  it. 

For  a  specified  history  of  loading  and  for  known  conditions  ahead  of  the 
front,  formally,  we  have 


eKL  eKl/XM’  5  eKL  eKL^XM’  ^ 

SKL  =  SKL^M5  ^  5  S  =  S^XM’  ’ 


(6.1) 


these  quantities  being  known  functions  of  position  and  time.  From  (6.1), 
we  can  compute  the  strain  rate  ahead  of  the  wave  front,  i.e., 

eKL  =  eKt/XM’  ^ 

Using  (6.1)  in  the  left-hand  side  of  (3.2),  we  may  then  define 

*(XM,  t)  =  t),  ej^,  t),  sCXjj,  t))  -  K  ,  (6.3) 

where  k  is  a  known  function  of  X^,  t  .  The  velocity  V  of  the 
elastic-plastic  interface,  and  hence  of  the  surface  of  acceleration 
discontinuity,  may  be  computed  from  (6.3)  by  a  formula  of  the  type 

(4.4). 

Since  e  and  V  are  determined  as  functions  of  X.  and  t  , 
lu-i  .  A 

equations  (4.23)  or  (4.25)  each  constitute  a  set  of  three  linear 
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inhomogeneous  equations  in  the  discontinuity  amplitudes  .  It  follows 
that  Xi  can  be  determined  uniquely  for  both  cases,  except  when  the 
determinants  of  the  coefficients  of  \i  vanish.  However,  these  cor¬ 
respond  to  cases  in  which  the  speeds  are  precisely  V/  ,  and  V/  \ 

(e)  (P) 

characterized  by  (4.16)  and  (4.19),  respectively.  For  example,  suppose 
the  elastic-plastic  interface  advances  into  a  plastic  region  with  the 
speed  satisfying  equation  (4.16).  Then,  provided  G.  4  0  ,  it 

follows  from  (4.23)  that  e^  =  0  .  Hence,  using  (3-30), 


•£*° 


(6.4) 


so  that  we  no  longer  have  a  region  of  loading  just  ahead  of  the  wave  — 
it  can  be  one  of  neutral  or  unloading  region. 

On  the  other  hand,  suppose  the  elastic-plastic  interface  advances 
into  a  region  of  unloading  or  neutral  loading  with  the  speed 
satisfying  equation  (4.20).  Then,  provided  0  ,  it  follows  from 

(4.25)  that 


(6.5) 


This  corresponds  to  propagation  into  a  region  of  neutral  loading,  if  we 
recall  (3-30).  Some  of  the  results  of  W.  A.  Green  [12]  obtained  under 
more  restrictive  conditions  are  compatible  with  those  given  here. 


drdsr  -Xves 


2.5  a  r^J-. 


?.  ^ui ss  or  &2 2  orders  riigr.er  ‘tns.n 


elastic  material  as  do  aoceLsratics.  distant inuit i e sT  He  shall  new  show 
that  identical  results  hmii  for  plastic  waves.  Since  by  (3-25)  and 
(3-23),  w,.  and  E^_.  are  functions  of  e,^,  e"_  and  S  ,  it 

fellows  free  tine  differentiation  of  (3.2?)  that 


“  _  Z  cZ  g  .  iff  o  -  o  Ov*  «  --  N 

bHL  ?o  =>3r-CiB  cs,_  ~  “AS  la,-,  '  cS 3  ^*KLMX“  "  UKL  'W 


•  PrAit?r»g  ~  ‘“21  rjS^As,v  "-,3  '  **a.K  5 


(7.1) 


during  leading.  Similarly,  free:  (3*29); 


-  f  -  C  .  C  .  o  C  1  A 

si*StCAS  ce._  ‘  A3  B<_  oS]  TOO 
A3  As 


i  /.  £  /  v  v  _  y  V  ^ 

‘  :'o  “SBBr  n»M  '  n,X  *  m,NJ  ’ 


(7*2) 


during  unloading  or  neutral  loading  - 

Consider  now  a  singular  surface  l(t)  of  order  3  with  re specs  to 
displacement!  Since  xm  x^,  e.^.,  e£_,  S  and  their  first  partial 


‘  See  -iruesdeli  [3]  for  a  general  derivation. 

r’fPor  a  singular  surface  of  order  it,  we  need  to  make  the  additional 
assumption  that  the  (m-2)th  time  derivative  of  F.  is  a  continuous 
function  of  X,  for  all  t. 


'  *,/  *  -I 
/  :  Y  ‘  <$£!/  ^ 


derivatives  are  continuous  at  such  a  surface  9  the  jump  of  the  time 
derivative  of  the  equations  of  motion  (2.6)  takes  the  form 

Xi,A  1  *A5,3  1  +  I  h,A3  1  *«  *  Ool^l  • 


(7.3) 


cl  ^ne. 
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side  of  s(t)  is  in  a  state  of  loading,  then  the  j ump  oi 
across  £{u)  is 

E  ®KL  I  =  po^lCM  “  C'KL  I  Xm,U  1  5 


(7.4) 


by  (?.l).  This  corresponds  to  the  case  (ii)  of  section  4  and  the  results 
corresponding  to  (i)  and  (iii)  can  be  discussed  in  a  similar  manner. 
Making  use  of  (4.7),  we  find  the  jump  relations  for  3rd  order  partial 
derivatives  of  displacement  and  2nd  order  partial  derivatives  of  stress . 
They  are 
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(7.5) 


(7.6) 


*>■’  ,  ** 


33. 


Here,  we  concern  ourselves  with  only  one  aspect  of  the  robL'-m,  namely 
the  determination  of  velocities  of  propagating  waves  ir  i  n  isotropic 
material,  using  the  results  of  section  4.  Thus,  we  ass  .je  that  a  state 
of  plastic  deformation  corresponding  to  (8.2)  and  (8.3)  .s  compatible 
with  the  field  equations  which  must  also  be  used  in  a  >  cnpletc  analysis 
of  the  problem. 

Before  proceeding  further,  we  need  to  recall  cert  .  n  results  con¬ 
cerning  the  forms  of  the  constitutive  functions  of  section.  3.  For  an 
initially  isotropic  material,  the  internal  energy  func  ;Lon  U  ,  the 
loading  function  f,  and  the  work- hardening  ten?,  r  h„  are 

ivij  ivLi 
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ftrcticz.  of  the  entropy  S 


anc  tne  -sen  -omt  invariants 
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(8.5) 


>*- —  :s^ly.  f  rust  be  expressible  as  a  function  of  S  and  the  ten 


cir.t  invariants  consisting  of  i",  if,  x"  and 

o  x5  2"  3 


-~L  SKK  5  "  S2CL  SKL  5  'S  SKL  SIM  SJsX  5 


(8.6) 


1.  =  e"  s_„  ,  L0  =  e"  e"  s._.  , 

i  KL  IjK  2  KL  IM  MK  5 


(3.7) 


“3  '  eKL  SLM  SMX  5  eKL  eIK  BMN  SKK 


Li.  =  e"  e"  s,_T  s. 


The  constitutive  eouation  for  8,,T  has  the  form  [l] 

AL 


Our  discussion  here  parallels  that  in  section  7  of  [l],  where  further 
details  can  be  found. 
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<t>  lew 


h12  “  h23  =  n3I  ~  °  5 


ar.d  the  components  of 


bf 
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KL 
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bf  o 
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^s12 

asa3 
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(8.14) 


From  the  results  (8.9)  to  (8.14),  we  can  now  obtain  the  following  informa¬ 
tion  about  the  coefficients  on  the  left-hand  sides  of  (3*25)  and  (3.28),  i.e., 
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*1112  "  A1113  "  A1123  "  A1213  “  A1222  “  A1223  *  A1233 


A1322  "  A1323  “  A1333  “  ^223  “  ^333  “  °  * 


C12  "  C23  “  C31  "  °  ’  D12  “  D23  “  D31  "  °  ’ 


H12  “  H23  "  H31  "  0  »  G12  "  °23  "  G31  "  °  * 


and 


A1212  “  A1313  “  ^323  * 

A1122  "  A1133  * 

*2222  “  A2233  “  A3333  ’ 

C22  =  C33  »  D22  *  D33  » 

**22  "  H33  ’  G22  “  G33  » 

B1212  =  B1313  “  B2323  * 

B1122  “  B1133  * 

B2211  =  B3311  * 

B2222  "  b2233  ■  B3322  “  B3333  ’ 


(8.15) 


(8.16) 
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together  with  results  analogous  to  the  first  of  (3.15)  for  B,„,nr  . 

With  the  help  of  (8.1),  (8.15)  and  (8.16),  the  value  of  and  the 

non-vanishing  components  of  ana  Hj  in  (4.17)2  and  (4.21)  are 

s(hT)  =  S11  A*l  +  s22^2  +  ^  » 

A11  =  Y  *-^1  Allll  +  ^  +  » 

^12  ~  ^  ^2  ^1122  +  ^1212^  5 

A13  ‘  *  h  K3(A1122  *  W  >  <8'17> 

A22  =  ^1  A1212  +  ^  +  W3^A2222  ’ 

^*23  =  N2  N3^A2222  +  A1212^  5 
A33  ^  +  K2^A1212  +  h2  A3333  5 

and 


G,  = 


V  N,  G. 


1  "11 


G«  = 


N2  G22 


G„  = 


n3g22 


(8.18) 


H,  = 


V  N,  H. 


1  “11 


*2  = 


N2  h,. 


H-  »  N. 


3  **22 


For  a  prescribed  wave  front  which  is  specified  by  its  normal  Nv 

K 

(generally  a  function  of  position),  we  may  substitute  (8.17)  and  (8.l8) 
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into  the  eigenvalue  equation  (4.20)  and  thereby  determine  the  wave 
speeds . 

Consider  now  a  plane  wave  front  whose  normal  is  at  an  arbitrary 
(oblique)  angle  to  the  X^-axis .  Without  loss  in  generality,  we  may 
choose  the  direction  of  the  X^-axis  to  be  perpendicular  to  the  plane 
of  the  X^-axis  and  the  normal  to  the  wave  front.  In  this  way,  the 
problem  is  reduced  to  a  two  dimensional  one,  so  that 

NK"(N1»N2»°)  »  (8-19) 

where  N^,  Kg  are  constants  (for  plane  waves)  satisfying 

+  K^  =  1  .  (8.20) 

In  view  of  (8.19),  (8.17)  and  (8.18)  simplify  and  the  non-vanishing 
components  of  (4.19)  become 

(Au  -  G1  H1  -  Y)X]_  +  (A1?  -  Gx  H2)\2  =  0  , 

(A12  -  G2  H1)X1  +  (A22  -  G2  H2  -  Y)X2  =  0  ,  (8.21) 

(a22  -  y)x3  =  0  , 

where  the  coefficients  in  (8.21)  are  obtained  from  those  in  (8.17)  and 
(8.18)  after  setting  <=  0  and  where  Y  is  defined  by  (5.10). 
Corresponding  to  the  eigenvalue 
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can  be  discussed  similarly  wi th  results  corresponding  to  (8.27)  and 
(5.2c)  in  the  iom 
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Xp  arbitrary 


Y  =  A11  ,  •  X-,  arbitrary 


(b  «30) 


Y  =  A2p  ,  X3  arbitrary 
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